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By a combination of analytical and numerical techniques, we demonstrate the existence of spherical, 
asymptotically flat traversable wormholes supported by exotic matter whose stress tensor relative 
to the orthonormal frame of Killing observers takes the form of a perfect fluid possessing anisotropic 
pressures and subject to linear equations of state: r = Apc^, P = fipc^ . We show that there exists 
a four parameter family of asymptotically flat spherical wormholes parametrized by the area of the 
throat ^(0), the gradient A(0) of the red shift factor evaluated at the throat as well as the values 
of (A,/i). The latter are subject to restrictions: A > 1 and 2/i > A or A < and 2^ < — |A|. 
For particular values of (A,/i), the stress tensor may be interpreted as representing a phantom 

' configuration, while for other values represents exotic matter. All solutions have the property that 

, the two asymptotically flat ends posses finite ADM mass. 
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I. INTRODUCTION 



The observed late time accelerated expansion of the universe [H [3] appears to imply that the universe is dominated by 
dark energy a substance having the property that if (pc^, P) stand for the energy density and isotropic pressure as 
measured by observes co- moving with the expansion, then P is negative or more precisely pc^ + 3P < 0. This property 
r , is often stated that the equation of state (EOS) for dark energy is w — < Additional observations of CMB, 

I^J^l gravitational lensing and supernovae mildly favor a more exotic EOS: w — < — 1 and such substance is dubbed 
as phantom energy. The possibility that the cosmic fluid is in state of phantom energy has another consequence. The 
stress tensor describing phantom energy violates the null energy condition and thus phantom energy is the type matter 
^ ' required to support wormholes [5;]. 

&^ ' The issue whether a static spherical distribution of phantom energy can support spherical wormholes has been 
0^ I addressed in [1, 0] ■ In these works phantom energy is modeled by a perfect fluid stress tensor possessing anisotropic 
' pressures with the radial pressure P,. = —t and energy density pc? obeying: w — = —-^ < — 1 while the 
. tangential pressure P is defined by the field equations. In Q an a-priori choice for the energy density has been made 
' whereas in 0] the wormhole metric has been a priori specified. 

, In the present work, we model exotic matter by a stress tensor so that relative to orthonormal frame of Killing 
observers takes the form of a perfect fluid possessing anisotropic pressures but we assume that pc^, t and P are 
subject to constraints referred as the EOS. Guided by the cosmological phantom EOS and generality purposes we 
adopt [i]: 



d' T = Xpc^, P — fipc^, X,p e Ft, (1) 

where (A, p) arc treated as free parameters subject only to the restriction that the stress tensor should violate the 
null energy condition in an open vicinity of the throat and restrictions arising by demanding asymptotic flatness. The 
problem of the existence of spherical wormholes supported by this type of exotic matter, is formulated as an initial 
value problem (IVP here after) with the throat serving as initial value surface (for more details regarding this approach 
see 0, [lo|). By a combination of analytical and numerical techniques we show that there exist initial conditions and 
values of (A, /z) where solutions of this IVP describe asymptotically flat non singular wormholes. Asymptotic flatness 
requires (A, /i) to obey either A > 1 and 2/i > A or A < and 2p < — |A|. Thus even though (A, /i) play a secondary 
role in the behavior of the solutions near the throat, they determine the asymptotic behavior of the solution. 

The structure of the present paper is as follows: In the following section we formulate the relevant IVP and discuss 
some properties of this IVP while in sections IlIII we discuss numerical solutions. In section HVl we present an analysis 
that explains the numerical outputs and advance arguments that set restrictions upon (A, p) so that the solutions 
describe asymptotically flat wormholes. We finish the paper with a discussion concerning the results and future work. 
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II. LOCAL WORMHOLES SUPPORTED BY A "PERFECT FLUID" 

We consider a static, spherically symmetric wormhole: 

g= -e2*Wdi2^d/2^r2(/)dn^ I e (-a, a), a > 0, (2) 

where r{l = 0) 7^ 0, ^^^^|o — 0, Jp'' |o > and thus I = marks the location of the throat. We support this wormhole 
with a stress tensor T^^ so that relative to the orthonormal frame of Killing observers decomposes according to: 

T<,p = p{l)c^u^ui3-T{l)XaXp + P{l){Y^Yp + ZaZp), (3) 



dt' dV r{l)de' r{l) sine df 

In the gauge of 0, if we define f{l) — r~^{l), introduce K{1) — f^^^^p' the trace of the extrinsic curvature of the 
50(3) orbits as embedded in t =cons hyper-surfaces and A{1) — then Gap = kT^p taking into account ([J) can 



be cast in the form : 



'^-^^^\K\l) + f^{l)~kp[l)c\ (5) 

- -mm) m + (1 - A + 2m)^, (6) 



dl 
dp{l) 



{\-im)-{i + ^)K{i) 



(7) 



-f^l) + + kp{iyx + A{l)K{l) = 0. (8) 

A throat of area A{0) ~ 47rr^(0) requires tension t(0) so that fcr(O) — r^^(O) and moreover: (t(0) — p(O)c^) > 
[Tl| . These conditions in view of ([T]) take the form: 

kXc'piO) - r-2(0), [r(0) ~ p{0)c'] = (A - l)p(0)c2 > 0, (9) 

and thus they hold, provided either (p(0) > 0, A > 1) or {p{Q) < 0,A < 0). If we adopt as a definition that an 
inhomogeneous phantom configuration satisfies w = — < ~1, then the first choice yields: w = — A < — 1 

and thus Tap describes a phantom like configuration at least in an open vicinity of the throat. For the second choice, 
even though p(0) < nevertheless t(0) is positive and for this case we interpret Tap as describing exotic matter jl^ . 

Any solution of describes a local wormhole having a throat of area A{Q) = 47rr^(0) = 47rf~^(0) at a prescribed 
A(0), provided it satisfies the initial conditions (Tlj |: 

/ 47r 47r 
V^(0)/ \kc^A{0) 

The system (|4][7]) combined with these initial conditions constitutes a well defined IVP. The theorem of "Picard- 
Lindelof" 13] assures the local existence of a unique solution defined on [~b,b] C (—a, a). Ideally, we would like 
to find necessary or (and) sufficient conditions so that these local solutions are extendible for all I E (—00,00) and 
moreover represent asymptotically flat wormholes. Due to the non linearities in (|3][7]) this is a formidable task. In this 
work we shall gain insights into the properties of the maximal solutions by first resorting to numerical techniques. It 
is our hope that these insights will eventually contribute towards an analytical treatment of the problem. 
For the purpose of numerical integrations, we compactify the coordinate I via: 

l(x) = -^, xei-1,1), ro = r(0). (11) 
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This transformation maps I ±00 into x — > ± 1, Z = into x = and transforms ([4][7]) into: 



{1-x 



dx 
2\2dK{x) 



-^-^{l + x^)K{x)f{x), 



dx 



(l^x-) 



2 2rfA(x) 



dx 



^^dp{x)c^ 



dx 



ro(l + a;2) 
ro(l+a;2) 
ro(l+a;2) 



--K'^{x)+f'^{x)-kp{x)c^ 



-K{x)K{x) - A2(x) + M^l£!(i - A + 2/i) 



i - 1] A(x) - (1 + ^) K{x) 



p{x)c^ 



(12) 
(13) 

(14) 

(15) 



and leaves the initial conditions (fTU]) form invariant. We integrate the IVP defined by these equations combined with 
(|10p for a variety of initial conditions and values of (A, /i). Before we enter into numerics, we mention three properties 
of smooth solutions of this IVP. 

• Any solution (f (x), iir(a:), A(x), p(x)) generated by initial conditions so that A(0) — 0, is "reflectionally 
symmetric" relative to the throat i.e. under " parity" transformation x ~x behaves according to p^ : 

r{x) = r{-~x), K{x) = -K{-x), k{x) = -A(-x), p(x) = p{~x). (16) 

• If {f{x), K{x), A(a;), p{x)) is any solution generated by an arbitrary set of initial conditions, then the functions: 



R{x) = r(-a;), Ki{x) -~K{-x), L{x) = -A(-a;), p{x) = p{-x), 



(17) 



define a new solution that satisfies the same initial conditions as {f{x),K{x),A{x),p{x)) does, except that 
L(0) = -A(0). 



If {f{x),K(x),A{x),p{x)) is any solution, then under rescaling of the throat 



it follows that 



f(0) f (0) = Af{0), A>0, 



Af{x), AK{x), AA{x), A'^p{x), 



(18) 



(19) 



satisfy (|12m5p and the rescaled initial conditions ^TU\i . 



Property (|16p implies that for reflectionally symmetric solutions, integrating (|12j|15p on the domain (0, 1) would be 
sufficient. Property ([TT]) allows us to concentrate on solutions obeying A(0) > while permits us to set the 
throat radius at some convenient value. 

III. NUMERICAL RESULTS 

We integrate I|12m5p subject to PU)) . using a Runge-Kutta integrator of different accuracy and in order to avoid the 
singularity at x = ±1, we stagger the exact points x — ±1. In all runs we employ units so that c = G = 1 and thus 
k = Stt. We set the throat radius r(0) = 1 and employ the following two sets of initial conditions (TTj : 



f(0) = l, X(0) = 0, A(0), p(0) = ^, A>1, 
f(0) = l, K{0) = 0, A(0), p(0)=g^, A<0. 



(20) 
(21) 
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At first we consider reflectionally symmetric solutions i.e. set A(0) = and choose A = 1.5 and /i G {1, 0.5, —0.5}. 
The resulting solution curves are shown in Fig. [T] In a second run we maintain A(0) = but choose (A ~ — 0.5,/i = 
0.5), (A = — 1, /i = —0.25) and (A = —1, fi = 0.5) and the graphs are shown in Fig. [21 

A number of runs aims to get insights in the behavior of reflectionally symmetric solutions upon changing the value 
of A(0). For these runs we consider a reflectionally symmetric solution for fixed (A, /i), and vary A(0) from A(0) — 
towards positive and negative values. We have used three different values for (A,/i) and the results of these runs are 
shown in Figs. I3I4I5I 



Case 


A 




Behavior 


1 


1.5 


1 


Asym. Flat 


2 


1.5 


0.5 


Decaying Non - Asym. Flat 


3 


1.5 


-0.5 


Non - Asym.Flat 


4 


-0.5 


0.5 


Non - Asym.Flat 


5 


-1 


-0.25 


Decaying Non - Asym.Flat 


6 


-1 


-0.5 


Asym.Flat 



TABLE L Table showing the values of the parameters (A, fi) used in Figs. I1I2I 



A(0)=0 A(0)=0 




-0.8 - 

-1 I 1 1 1 ' 1 I 1 1 1 

-1 -0.5 0.5 1 -1 -0.5 0.5 1 

X X 

FIG. 1: Reflectionally symmetric solutions generated choosing: A(0) = and (A,^t) labeled as case (1,2,3) in Table I.) 

IV. EXISTENCE OF ASYMPTOTICALLY FLAT WORMHOLES 

The numerical outputs displayed in Figs. [l][S]show that solutions are divided into two families: the first one contains 
solutions where all variables decay to zero as x ^ ± 1 whereas the second one contains solutions where the variables 
become unbounded in one (or both) ends. This behavior raises a number of questions: Which if any, of the decaying 
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A(0)=0 



A(0)=0 



< 









X X 

FIG. 2: Reflectionally symmetric solutions generated by choosing: A(0) — 0, and (A, /i) labeled as case (4,5,6) in Table.) 

to zero solutions describe asymptotically flat wormholes? How the initial conditions and values of (A, /i) affect the 
global behavior of the solutions? Bellow, we provide an explanation of that behavior and draw a few conclusions 
based on these graphs. 

We begin by first discussing an exact solution of ([4][8l). This solutions is generated by choosing (A,//) so that 
1 — A + 2/Lt = and taking A{1) = as the global solution of eq. The restriction 1 — \ + 2fi — implies that the 
EOS takes the form: 



T = \pc^, P = ^ ^ ^ pc^, A > 1 or A < 0. 
and the resulting solution in curvature coordinates is described by (for details see (isj'l: 



(22) 



kp{r)c^ 



^2A V r 



3A-1 



— + r dQ , r £ (ro,oo), 



A 



, r e (ro,cx)). 



(23) 
(24) 



For latter use, we notice that for A > 1, the density decays according to kp{r)c? = 0{r '^+'^-'), < e < 1, while for 
A < 0, kpc^ = 0(r-(3+e))^ g > Q 

In the remaining part of this section we analyze the asymptotic behavior of the solutions. We consider a solution that 
decays to zero as / — s- ±oo (the graphs in Figs. [T][5]show that such solutions exist). For such solution and away from 
the throat we introduce coordinates {t, r, 9, (p) so that: 
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^2, |a,=0.5 ^2, |a,=0.5 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 




FIG. 3: In these graphs we have taken A = 2, /j. — 0.5 so that 1 — X + 2fi — 0, and we have kept them constants throughout. 
We graduaUy change A(0) starting from A(0) — and increasing (decreasing) towards positive (negative) values. Notice that 
after a characteristic value of A(0) the solutions become unbounded at one end. 



g = -e^^Wdt^ + + rW, r>ro= r(0) > 0, (25) 

r 

m{r) = ? + ^ f r^p{r')dr'. (26) 



ro 



2 2 

Either from (|3][7]) or directly from Gap — kTap we find 



= — ^ pW, r e (ro,oo), (27) 
dr{r) ^/^2 , m d<i>{r) 2[P(r) + r(r)] 

d$(r) -krMr^ + 2m(r) , , , , 

^ - ^ ' ^ \ r e (ro,«3). (29) 



dr 2r{r — 2m(r)) 

Taking into account the EOS and introducing the variables: 



W{r) = R{r) = kpiryr^, (30) 



(l?7ll^ become: 
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h=l.5,\x=l 



^1.5,|a=l 




A(0)=0.8 
A(0)=0.6 
A(0)=0.4 ■ 
A(0)=0.2 - 
A(0)=0.0 ■ 
A(0)=-0.2 
A(0)=-0.4 ■ 
A(0)=-0.6 
A(0)=-0.8 ■ 












A(0)=0.8 




A(0)=0.6 


#■•/ \ 


A(0)=0.4 


#7 \ 


A(0)=0.2 




A(0)=0.0 




A(0)=-0.2 




A(0)=-0.4 




A(0)=-0.6 




A(0)=-0.8 

















'h=l.5,\i=l 



'h=l.5,\i=l 



A(0)=0.8 
A 0)=0.6 
A 0)=0.4 ■ 
A 0)=0.2 ■ 
A 0)=0 ■ 
A(0)=-0.2 
A =-0.4 • 
AO =-0.6 
-A(0 =-0.8 





FIG. 4: In the graphs we have taken A = 1.5, ^ = 1. Like in Fig. ((3}, we vary A(0) from A(0) = up to A(0) = 0.8 and 
A(0) = —0.8. Again after some characteristic value of A(0) the solutions become diverging at one end. 



dW{r) 
dr 



-W{r)+R{r), (31) 



(33) 



dr \ ^ ' \\ / 2(1 - W{r)) 

d$(r) _ -\R{r) + W{r) 
dr ^ 2(1-M^(r)) ' 



This system becomes singular at three places: as r ^ tq, on any r S (rp, c») so that W{r) = 1, and as r ^ oo. Since 
for our background solution all fields decay to zero as I ±oo, we assume that pm33p is defined in a domain [i?o, oo) 
where i?o is sufhciently large so that: 

W{r) <l, R{r) <l ,\/ r e [i?o,oo). (34) 
and additionally: lim.r^oo W{r) = 0, limr^oo Ri''') = 0. On [i?o, oo) we introduce a new variable t via: 

r = roe*,t € [T,oo), T = log f (35) 



so that pHI33p becomes an autonomous system: 
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dW{t) 
dt 



-W{t) + R{t), (36) 



(38) 



dt X ^ ' \X J 2{l ~W{t)) 

d$(t) _ -~XR{t) + W{t) 
dt ^ 2{1 -W{t)) 



Since solutions of ((38)) are determined by {W{t), R{t)), we restrict our attention to the first two equations. They can 
be written in the form: 



dt 

where x{t) = {W{t), R^t j)-^ , A stands for the matrix 



"^""^^^ - Ax(t)+r(x(t)), t e [T,oo), (39) 



while F is defined by: 



F:AclR^^lR^:iW,R)^F{W,R)^0, (^-1) ^^ ^^ fl. (41) 
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According to (p4)) . we take: 

A={^={R,W) I |x|<i?<l}, (42) 

where we employ: |x| = + |a;^| as the norm of M^, JJj. On this A, it is easily seen that F is actually analytic, 
and moreover it satisfies: 

|F(x)| < fM^i-^ Ixl) |x| < iV|x|, (43) 



where 



|A| 



fl — A)^ 

M = max weA \[2{1-W)r^\ and N = max ^ g a ( M , |x| 



Let now the IVP: 

dx(t) 



|A| 



= Ax(i) + r(x(t)), t G [T,oo), (44) 



dt 

x(T) = xo, xo e A. (45) 

By the Picard-Lindelof's theorem [isj . there exist unique solutions x(t) defined on [T, e), e > 0, and by the variation 
of constants formula this solution is described by: 

x(i) = e(*-^)^xo + /*e(*-^)^r(x(s))ds, t G [T,e). (46) 

Although this x(<) is defined only on [T,e), e > 0, as long as |xo| is small enough and under some conditions upon 
the eigenvalue Ai — —1, A2 = — ^ of A, it can be continued so that remains in A V t G [T, 00). For that, let us 
suppose (A,/i) are chosen so that Au > 0. This restriction implies A2 = — ^ < and thus both eigenvalues of A are 
negative. Standard estimates [H, [i3| show that there exist K > and p > max(Ai, A2) such that: 

|e(*-T)A| < ^gP(t-T) V t > r. (47) 

As long as max(Ai, A2) < 0, we can always choose a p so that: max(Ai, A2) < — |p| < 0. For such p we obtain from 
dlSl) and 621): 

|x(t)| <ife-l''l(*-^'|xo|+^y" e-l''l(*-")|F(i(s))|ds. (48) 
Upon multiplying both sides by e''''* we arrive at: 



e\p\*\^{t)\<Ke\p\^\^o\+KNf el''l''|x(s)|ds, (49) 

Jt 

where we made use of (l43l) . By appealing to Gronwall inequality [3] we conclude that: 

|x(t)| < i^lxole-^*-^)!!"!-^^!, t > T, (50) 

and by shrinking the size of A if necessary, we can always make \p\ — KN > and thus the solution x(t) of (|44l45p 
remains in A, V t > T. 
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This analysis shows that as long as A/i > 0, the solutions to the IVP (|44l45p are dominated by the linear part and 
since both eigenvalues (Ai, A2) are negative, they decay exponentially to zero as i — *■ 00. To get more insights on their 
behavior we consider: 



W{r) = - + -^, e>0, (51) 



r 



and returning to pOII33p we find: 



W{r) = - + (52) 



Pir) = + o(^), Po^ -(^ - l)ci, (53) 



^ir) = 1>o - f - £l^(i±^^ 1 + 0(1). (54) 



2r 4/i r'^'x r 

Always under the assumption A/i > 0, if A2 < Ai = — 1 then ^ = l + e, e>0 and thus 



W{r) = ^ + ^, (55) 

Pir) = -^^i^ = -;S7>^>0, (56) 
c ciA(l + 2fi- X) 1 



If Ai < A2 < then ^ = 1 — e, e>0 and the solutions exhibit slower decay rates: 



W{r) = - + e > 0, (58) 



p{r) = = e>0, (59) 



c ci(l + 2/i-A) 1 
2r 4/7 



It is interesting to note that for the particular case where 1 — A + 2/i = 0, we obtain A2 = — ^ = ^ 1 + j and thus 
if A > 1 then Ai < A2 < 0, which implies kpc^ = 0(r"(^+'')), while for A < A2 < Ai, and thus kpc^ = 0{r~'-^+^^). 
This behavior is in agreement with the exact solution described by (|23I24[) . 

There is an important difference between the asymptotic behavior of solutions described by ()55II57|) and those 
described by (|58ll60p . The first family has finite ADM mass Madm while for the second family this mass is actually 
divergent. In order to see that, we notice that the spatial metric ^'^^7 can be written in the form: 

^'^^ = T31F(;:) +^'^"'' ^ (^0,00), (61) 

and by employing conformal coordinates takes the form: 

(3)7 = 9?{R)[dR^ + R^dn^] = 9?{x)[dx^ + dy^ + dz\ (62) 

where 



~2 (dr{R)V 1 



In order to evaluate the Madm, we employ the representation 
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FIG. 6: Sketch of the parameter space (A,^). The regions labeled by (3,4) do not generate solutions that decay to zero 
asymptotically. Regions labeled (1, 6) generate solutions decaying to zero at infinity but the decay is slow. Finally regions (2, 5) 
generate solutions that decay sufficiently fast so that they have finite ADM mass. 



Madm = / ['^(7) + 77'""7"''7'''(27mn.<i7ac,b - 27,„„,rf7„c,6 (64) 

167r 4 

+7dm,67nc,a - ^cd,almn.b)\^/l<f' X , 

where E stands for any asymptotic end defined by restricting r so that: r > Rq. 
By appealing to the Hamiltonian constraint ^^"^ R{'-f) — 2kpc^ and (l62l) . the right hand side of (|64p yields: 



Madm = 



2kpc^ - 



2tl^{R) \ dR 



(65) 



However it can be easily seen that this integral converges for the solutions described by (|55|56p and diverges for those 
described by (|58l59p (for more details see also [§]). 



V. DISCUSSION 



The results of the present paper can be succinctly summarized by partitioning the (A, p) plane according to whether 
the eigenvalues (Ai,A2) of the matrix A satisfy: A2 < Ai or Ai < A2 < (see Fig. |6]). The values of (A,/x) required 
to generate asymptotically flat solutions in the sense that the Madm is finite is indicated. Our analysis shows that 
solutions of (|12j|15p subject to (|20l2ip if they are decaying as r —> cx3, then (A, /i) should lie in the regions (1, 2) and 
(5,6) of Fig. inland the decay rates are those described in (|55ll57p and (|58ll60p . It is important however to stress 
a point. Although the numerical outputs show that there exist initial conditions so that the solutions are decaying 
as r — > 00 and moreover are asymptotically flat whenever (A, /i) are taken in regions (2,5), we have not shown that 
for any (A,/i), lying on regions (2) or (5) of Fig. [SI there exist initial conditions so that the solutions of (|12m5p are 
reaching the asymptotic region. This for the moment is an open issue. 

Our results establish the existence of reflectionally symmetric asymptotically flat wormholes having throat of arbi- 
trary area. This conclusion is a consequence of the rescaling property of the solutions under rescaling of the throat 
area. As far as non reflectional wormholes are concerned, the graphs in Figs. |3][n] show that they do exist (notice 
for these figures A2 < —1). However they do not exist for all values of A(0). The numerical outputs show that after 
some critical values of A(0) the solutions become unbounded. This suggests that A{0) and A(0) may not be taken 
independently although the nature of any constraint of this sort is for the moment unknown [20l |. 
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